Properties of two - dimensional dusty plasma clusters. 
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Two-dimensional classical cluster of particles interacting through a screened Coulomb potential 
is studied. This system can be used as a model for "dusty particles" in high-frequency discharge 
plasma. For systems consisting of N — 2 — 40 particles and confined by a harmonic potential we find 
ground-state configurations, eigenfrequencies and eigenvectors for the normal modes as a function 
of the Debye screening length Rd in plasma. Variations in Rr> cause changes in the ground-state 
structure of clusters, each structural rearrangement can be considered as a phase transition of first 
or second order (with respect to parameter Rd)- Monte Carlo and molecular dynamics are used 
to study in detail the melting of the clusters as the temperature is increased. By varying the 
density and the temperature of plasma, to which the particles are immersed, one can modulate 
thermodynamical properties of the system, transforming it in a controllable way to an ordered 
(crystal- like), orientationaly disordered or totally disordered (liquid-like) states. The possibility of 
dynamical coexistence phenomena in small clusters is discussed. 
PACS: 64.70.Dv; 73.20.Dx; 61.20.J 



I. INTRODUCTION 



Small charged particles of "dust" are rather common systems, and are observed on different scales and in different 
environments: clusters of dust in the interstellar medium, charged colloidal suspensions, ordered structures in the gas 
discharge used in thermal processing of materials are examples of such systems |^-|^] . At present much attention is 
paid to experimental investigation of the properties of " dusty plasma" which is a system of small micrometer particles 
in a high frequency gas discharge. One of the main reasons of this attention paid to the artificial objects like this one 
is the ability of a direct observation of their static and dynamical properties. The study of dusty plasma crystals and 
liquids 'in vitro' which is being carried out in a number of laboratories around the world Q is of a great importance for 
understanding the plasma properties and is a powerful tool for the examination of melting, annealing and formation 
of defects of different kinds. 

Small particles immersed in a plasma may acquire large (up to 10 5 e) charges Ze due to high mobility of plasma 
electrons. The presence of plasma screening modifies the Coulomb interparticle interaction and, with a good preci- 
sion Mffy , the system can be described as a system of particles interacting with a Yukava-type pair potential. Here we 
consider two-dimensional (2D) clusters of dusty particles confined by an external harmonic potential of strength a. 
It is obvious that it is 2D system that realizes when one consider particles of 'dust' immersed in a plasma discharge 
cloud with the transverse dimension higher than the Debye screening length. The energy of the system has the form: 

E^iZe^ ^tMM + a^ (1) 

Kj l r «l i=l 

— 1/2 

The Debye screening length in plasma is determined as Rrj = [Anqfrii/kbTi + 4ire 2 n e /ki,T e ) where (foTij and Tj 
are the charge, mean density and temperature of plasma ions and e, n e , T e are that of plasma electrons respectively. 

The energy of a cluster, being written in dimensionless units ro = (Ze) 2 / 3 /a 1 / 3 for distances and Eq = ar\ for 
energies becomes: 

E = £ !2 P(z7l* i D + £ |ri|2 (2) 

i<j 1 lJ ' i=l 

where the dimensionless parameter 7 = rg /Rd defines the range of the pair interaction potential. From (^|) one can see 
that the thermodynamic state of a cluster of a given number of particles is determined by two dimensionless parameters: 
the inverse dimensionless screening length 7 and the dimensionless temperature of "dusty" grains = kbT ' j 'Eq. The 
range of interaction between particles in a cluster 1/7 = t/Rd is controlled by the density and the temperature of a 
plasma (see above). 

In this paper we consider the properties of two-dimensional (2D) clusters (|J) as a function of the number of particles 
N < 40, the screening length I/7 and the temperature 0. We show (see Sec. |0|) that the change in the screening 
length (i.e. in parameter 7) causes rearrangements of the ground-state structure of the cluster at a set of points 7*, 
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such structural transitions can be treated as phase transitions of the first or the second order with respect to parameter 
7. In Sec. II] we apply molecular dynamics (MD) and Monte Carlo (MC) simulations in a canonical ensemble in 
order to study the thermodynamic properties of small clusters. We show that in clusters of rather a small number 
of particles and at a small enough plasma screening (at 7 < 10), as the temperature is increased, the orientational 
disordering happens first, i.e shells rotate with respect to each other by losing their mutual orientation order. At more 
higher temperatures a total disordering of cluster shells takes place. 



II. GROUND - STATE CONFIGURATIONS 

Ground-state configurations (see Table 1 and Fig. 1-3) of the system (^) have been found with the help of the 
following methods: 1) The modified Newton method ||; 2) The combination of "random search" and "gradient 
search" methods For the results to be more reliable, all configurations discussed below have been independently 
obtained by both of this methods. Of course, no one of the present methods of the minimization of a multidimensional 
function is able to guarantee that the configuration obtained is a global minimum one. To overcome this difficulty we 
have used up to 200 randomly distributed initial configurations. This approach have also enabled us to investigate 
both the local minima and their caption regions (i.e. "specific weights" of local minima). 

At 7 <C 1 model (^) describes the Coulomb cluster in a harmonic trap, the system that have been actively studied 
both experimentally Jn],[llJ and with the use of computer simulations jl^j , j|,[l3|. In particular, the calculations 
carried out before have revealed that particles in small finite systems arrange themselves into shells. An analysis of 
shell structures for different number of particles N enables one to consider the system as belonging to some period 
of a Mendeleev - type table. This table can be viewed as a classical equivalent to the well - known Periodic Table of 
elements. 

The presence of parameter 7 which determines the range of the pair potential enables one to investigate the influence 
of this range on the structures and properties of clusters. The fact that cluster structure depends on the range of 
interaction potential becomes obvious from Table 1 in which some ground-state configurations for 2D clusters in a 
harmonic trap are presented. 

As the value of parameter 7 changes, rearrangements of ground-state structure take place, each point 7* of any 
of these chan ges can be treated as a point of phase transition of one kind or another. Following by the approach 
used in work pq[ , the order of these phase transitions can be determined from the plot of the ground-state energy 
-E(7), the discontinuity in the n — th derivative of E(j) with respect to parameter 7 corresponds to the phase 
transition of n — th order. Another way to determine the order of the phase transition is to analyse the behaviour of 
eigenfrequencies ^(7), i = 1, 27V for the normal modes: first order transition takes place at the point 7* at which any 
of the eigenfrequencies exhibits a jump while softening of any of the eigenfrequencies (when it becomes zero) testifies 
about second order phase transition. 

The eigenfrequencies of the cluster of N = 10 particles vs. screening parameter 7 £ [0, 10] are presented in Fig. 1. 
At points 7 « 1.4 and 7 ~ 8.2 of first order transitions the eigenfrequencies exibit jumps that are clearly seen. 
From Fig. lb one can see that with a decrease in the interaction range first, at 7 w 1.4, the distribution of particles 
throughout shells changes and the configuration typical for the Coulomb interaction is replaced by one appropriate to 
the dipole cluster of 10 particles ({2, 8} — > {3, 7}). Further reduction in the screening radius transforms (at 7 w 8.2) 
the cluster to the most "packed" state {2, 8} which is characteristic of a system of hard spheres. 

In Fig. 2 the ground-state energy of a cluster of 33 particles are given. At 7 « 3.751 the first derivative of the 
ground-state energy with respect to parameter 7 is discontinuous (see inset of Fig. 2a). Investigation of cluster 
configurations shows that the numbers of particles in two outer shells change here as {1, 6, 11, 15} — > {1, 6, 12, 14}. 

The point of the first order phase transition 7* can be determined as the point, where the energies of ground- 
and the lowest excited (metastable) states are equal. This statement is illustrated in Fig. 2b in which both ground- 
state energy E{pf) and energy E^(-j) of the lowest local minimum are plotted in the vicinity of the transition point 
7* w 3.751. From this figure one can see that the configuration of the global minimum at 7 < 7* corresponds to a 
local one at 7 > 7*. 

One can see that in the cluster of 37 dipole particles one of the particles is between the second and the third shell 
to form an interstitial (analogous to the Frenkel defect in crystals) and to make the division of the ground-state 
configuration into shells ambiguous ||. The 'dusty' cluster of 37 particles is supposed to exhibit a rich variety of 
structural rearrangements while varying 7. The investigation of this cluster at different values of screening strength 
has revealed four phase transitions in the region 7 £ [0,1.6] (see Fig. 3a), namely two second order transitions (at 
7 w 0.78 and 7 rs 1.22) and two first order transitions (at 7 R* 0.52 and 7 R* 1.34). 

From Fig. 3a on can see that the number of particles in the outer shells changes at 7 w 0.52. It is worth while 
to note that a distinctive feature of the first order phase transition is the abrupt change in the cluster structure, 
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Usually, this peculiarity exhibits as the change in the distribution of particles throughout shells (like that one can 
observe for the clusters of 10 and 33 particles, see Fig. 1,2). This very change takes place at 7 = 0.52 for the cluster 
involved. However no apparent changes in structure of the cluster are seen at the point 7 ~ 1.34 of the first order 
phase transition. More detailed study have shown that at this point there exist a rotation of the third shell with 
respect to the fourth one. This is illustrated in Fig. 3b which presents the mutual orientational order parameter g SlS2 
of different pairs of shells {si, S2} ||, the value which is very sensitive to the changes in mutual orientation of cluster 
shells. 

Subsequent increase in the parameter 7 leads to the realization of two other first order transitions which are depicted 
in Fig. 3b. In the first of them (at 7 7.015) one of the particles implants between the second and the third shells (see 
Table 1 and the discussion above). The corresponding transition can be written as {1, 7, 13, 16} — > {1, 6, 1, 12, 17}. 
At 7 > 19 the cluster becomes well - facetted and has the most symmetrical structure {1,6,12,18}. Note that in 
this region of 7 the minimal nonzero eigenfrequency Lu m i n corresponds to twofold degenerate vibrations of the whole 
cluster in the harmonic trap with the frequency oj m i n — \/2. Further decreasing in the range of the pair potential 
does not lead to any structural rearrangements. 

The study of Coulomb and dipole clusters have shown that the basis for most configurations is provided by different 
parts of 2D hexagonal lattice ]g|,13|]. When describing and analyzing the properties of such configurations it is suitable 
to introduce into consideration the "crystal shells" Cr c that are concentric groups of nodes of ideal 2D crystal with 
c nodes placed in the center of these groups. Obviously, in view of the axial symmetry of the confinement potential, 
we can concentrate on a finite number of the most symmetrical crystal shells which, by the number of particles in the 
center of the system, can be divided into the following groups: CVi, CY 2 , Cr 3 , CY4. With the help of the crystal shell 
concept we have found that changes in the ground state structure of "dusty" clusters, as parameter 7 is increased, 
comes in such a way, as to fill the maximal number of crystal shells. 



One of the distinctive peculiarities of small clusters is the existence of two stages of their disordering |l^,[l4j : an 
intershell (orientational melting of shells si and S2 at the temperature 6 SlS2 ) and a radial disordering (total melting 
at temperature 9/). The analysis of eigenfrequencies shows that the clusters with the small values of lowest nonzero 
eigenfrequencies uj m i„ have the eigenvectors corresponding to mutual rotations of cluster shells. Such clusters have low 
temperatures 6 SlS2 °f intershell disordering. It is obvious that changes in the cluster structure caused by variations 
in the control parameter 7 lead to the modulation of the temperatures of both orientational SlS2 and total 0/ 
disordering. Moreover, the phenomenon of the orientational disordering may disappear at all if the cluster has a well 
packed structure. The results of our simulations have proved this suggestion. 

The dependencies of the mutual orientational order parameter 321(6) for two-shell cluster of N — 10 particles 
at several values of parameter 7 are given in Fig. 4a. It is evident that g Sl s 2 drops to zero at the point of relative 
disordering (mutual rotation of shells s\ and S2) |l(|. One can see from Fig. 4a that the change in the system 
configuration {2,8} — > {3,7} which occurs at 7 » 1.4 (see Fig. 1) leads to a sharp decrease in the orientational 
disordering temperature: ©21(7 < 1.4) w 1.3 ■ 10~ 4 -> 621(7 > I- 4 ) ~ 0.7 • 10~ 5 . 

The cluster is well - packed in the region 7 > 8.2 (see Fig. 1) and that is why it does not experience the orientational 
melting, when an increase in the temperature leads directly to the interchange of particles between shells at 6 « 10 . 
This can be seen from the analysis of radial square deviations v%, : 



The dependence u^(6) is given in Fig. 4b, also shown are analogous curves at 7 = 1 and 7 = 2. One can see that 
even the slightest variation in the value of control parameter 7 may change the temperature of the total melting up 
to orders. 

The changes in the interaction potential lead to the modification of the structure of the energy surface which 
determines the type and the distinctive features of the phase transitions. For this reason, one can suppose that 
at some values of the parameter 7 the system can have very interesting thermodynamic properties. In Fig. 5a the 
dependence of the radial square deviations (||) of four - shell cluster of N = 33 particles at 7 = 3.76 is given. The 
graph has a number of plateaus located in different temperature intervals. A detailed investigation has shown that 
the regions of a sharp increase in u^. correspond to the radial disordering of different pairs of shells: particles start to 
interchange between the third and the fourth shell at temperature 6g 4 « 5 • 10~ 4 and between the second and the 
third - at O23 ~ 0.005. The total melting of the cluster takes place at 0? w 0.01. 



III. PHASE TRANSITIONS 
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Some useful information about the character of the disordering considered can be obtained by exploring the local 
minima distribution p(Ei oc ). In order to estimate this histogram, at each measurement time point we have performed 
several hundreds of gradient search iterations to find the nearest local minimum with the energy Ei oc . 

In Fig. 5b the local minima distribution of the system of 33 particles at 7 = 3.76, = 10~ 4 , 9 = 8- 10~ 3 is 
shown. In an entirely ordered state (at = 10 ) the system lives in the vicinity of the global minimum (with energy 
E = 64.795946 and the structure {1,6,12,14}). At = 8 • 10~ 3 the cluster can be also found at the lowest local 
minimum E^ = 64.795975 with the configuration {1, 6, 11, 15} (see Fig. 2b). 

Considering the results stated above one can conclude that the first disordering seen in the temperature interval 
G [10~ 4 , 10 -3 ] (see Fig. 2) corresponds to the nonzero probability of the cluster to be found in the state {1, 6, 11, 15} 
which is metastable at the given value of the parameter 7. Such changes in the distribution of particles throughout 
shells demand the overcoming of potential barrier that, knowing about large specific weights of both "ground" and 
"excited" states, allows one to treat this temperature interval as that of the dynamical coexistence of two cluster 
forms {1,6, 12,14} ^ {1,6, 11,15}. 



IV. CONCLUSIONS 



In this letter we have presented the results of a study of finite "dusty plasma" particle system. As a function 
of Dcbye screening length Rr> (for the particle charge in the plasma) we have found ground-state configurations of 
clusters consisting of N < 40 particles, their normal modes eigenvalues and the corresponding eigenvectors. The 
clusters undergo structural transitions which manifest itself as phase transitions of first or second orders with respect 
to parameter Rjj. At points of first order transitions cluster coordinates experience jumps which lead either to the 
change in the shell distribution or to the rotation of some shells relative to each other. At points of second order 
transitions one of the eigenfrequencies softens and particle coordinates change continuously. 

By varying Rr> (for example, by varying the temperature and the density of the plasma) one can modulate thermo- 
dynamic properties of the system and considerably change the temperatures of both the orientational and the total 
disordering. It have turned out that for some clusters, as the screening becomes sufficiently high, the disappearance 
of the orientational melting of different parts of the system takes place and an increase in the temperature leads 
staightway to the interchanges of particles between shells. 
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[16] the quantities tp s and g 31S2 = V'siV's-j are analogous to the parameter tpe and the correlation function g&(r) of infinite 2D 
systems, the softening of the correlation function ge(r) — > 0, r — > oo corresponds to the relative orientational disordering 
of remote parts of the system (with the absence of the translational order). 

Tabic 1. 

Ground-state shell structures {Ni, N 2 , ■■■} for dipole, Coulomb and logarithmical clusters of N particles confined in a 
harmonical potential. 



TV 


Dipole cluster 


Coulomb cluster 


Logarithmical cluster 


9 


2,7 


2,7 


1,8 


10 


3,7 


2,8 


2,8 


11 


3,8 


3,8 


3,8 


32 


1,6,12,13 


1,5,11,15 


4,11,17 


33 


1,6,12,14 


1,6,11,15 


5,11,17 


34 


1,6,12,15 


1,6,12,15 


1,5,11,17 


36 


1,6,12,17 


1,6,12,17 


1,6,12,17 


37 


1,6,1,13,16 


1,7,12,17 


1,6,12,18 


38 


2,8,13,15 


1,7,13,17 


1,6,12,19 



Fig. 1. 

Eigenfrequencies (a) and the lowest nonzero eigenfrequency Lu m in (b) of the 'dusty' cluster of 10 particles. Inset: 
ground state configurations in three different regions of control parameter 7. 



Fig. 2. 

Cluster of 33 particles, a) The first derivation of the cluster ground energy with respect to 7. In the inset the region 
of the first order phase transition is shown on an enlarge scale, b) Energies and configurations of the lowest local 
minimum (with the energy measured from the ground state energy E) of the cluster in the region of the phase 
transition. 



Fig. 3. 

Cluster of 37 particles. ) The lowest nonzero eigenfrequency 0J m in{l) and mutual orientational order parameter of 
different pairs of shells g Sl s 2 {l)- b) The regime of strong plasma screening. Eigenfrequency oj m in corresponds at 
7 > 19 to the motion of the cluster at a whole in the confinement. 



Fig. 4. 

Two - shell cluster of 10 particles. ) Thermodynamical mean of the mutual orientational order parameter < 1721 > (@) 
at different values of 7. b) radial mean - square deviations od particles vs. temperature uf(Q). 



Fig. 5. 

Four - shell cluster of 33 particles a) radial mean - square deviations, b) local minima distribution histogramm 
p{E^ loc ^) of the cluster in the ordered state (at 6 = 10~ 4 ) and at 9 = 8 • 10~ 3 , when there is an interchange of 
particles between the third and the fourth shells. 
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